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The Fast Fourier Transform (FFT) algorithm has 
keen the omcst attractive procedure tc carry cut 
digital Ciscrete Fourier Transform (DFT) analysis. 
However, interest has been generated in two ctker 
algorithrs, the chirp-z-transform (CZT) and the prime 
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develcred. Their results were compared with results 
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analysis tc parametric variations of their hardware 
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impulse, rectangular, and composite sinuscidal input 
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I. INTRODUCTION 


— ae ee ee Se ee ee ae 


A. EACKCECUNT 


The Fcurier transform has been with us many years and 
the uniqueness of the transfcerm allows Fourier analysis 
technigues developed in one field to be agfplied to many 
diverse areas. Typical application areas of the Fceurier 
transforma include: 


S -~ The “Fourier transfcrm of the 
near system 1S given by the croduct 
m transfer function and the Fourier 
he input signal. 


Antennas — The field fattern cf an antenna is 
given by the Fourier transform of the antenna 
current illusination. 


Qgtics -— Cptical systems have the peo heeey that a 
Fourler transform relation exists kLetween the 
ligkt anoplitude distribution at the frernt and Lack 
focal planes cf a converging lens. 


Randcm Ercecess — The power density spectrum of a 
randcm process 1s given by the Feurier transfcre 
cf tke atto-ccrrelation function cf the process. 


Probability -— The characteristic furction of a 
Candcm variable is defined as . the Fourier 
transform cf the probakility density function of 
the randc@ variable. 


Cuantum Fhysics — The “uncertainty principle" in 
guantus teeny is fundamentally asscciate with 
he. Fcurier transform since particle scmentum ana 
pesiticr are essentially related through the 
FCOurler transforg. 


Eounéary-Value Problems — The solution of partial 
differential equations can te obtained Ly means of 
the Fetrier transform. 





Althcugh Fourier analysis allows cne to more easily 
examine a functicn from another point of view,i.e., the 
transform dcmuain, the methcds available before 1665 to 
ccgpute tke Feurier ccefficients were very time consuming 
and ccstly. Then in 1965 Cocley and Tukey fublished their 
wathematical algorithn which computed the Fcourier 
coeificents with such less computation effort than haé keen 
required in the past. This methced has become known ag the 
"fast Fetrier transform" or FFT and has produced pwajor 
changes in ccmputational techniques used in digital spectral 
analysis, filter simulation and related fields. 


Eieneugh the FrYT is the most well knewn algoritbhie tc 
ccupute the Fcurier coefficients, it is net the only trethod 
fer computinc the "discrete Fourier transform" (DFT). The 
chirp-z-transfcorm (CZT) is an algorithm for evaluating the 
z-transfcrm cf a finite duration sequence along certain 
general ccntcurs in the z-plane. Evaluating the DFT cf the 
sequence is a special case of the CZT, i.e., the z-transform 
cf the finite duration seguence is evaluated on the unit 
circle in the z-plane. Although the C2ZT is not quite as 
efficient as the FFT, it eliminates many of the restrictions 
of the latter. These restrictions are listed at the end of 
the CZT alccrithm derivation. The “prime transform" is an 
Sroecritumecndt fcllows directly frem the CZT derivatior. ft 
is based ufen "N", the number cf samples cf data, rEeing an 


cdd pride. 


The CZI and prime transform algcrithms are receiving 
considerarle attenticn for applications in spectral 
analysis. this has been brought about because the bulk of 
the computaticns of these algorithms is ferformed Ly a 
transversal filter, which is easily isplemented Ly charge 
transfer devices (CID) or surface acoustic wave (SAW) 


devices. Haxréware and software implementaticns of the FFT, 
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CZT, and fripwe transform algcrithms are discussed later in 
this charter. 


Be. EASIC CCNHCEFTS OF FCURIER TRANSFORM ANALYSIS 


The esserce of the Fourier transfcrm of a waveform is to 
deccmposeé or separate the waveform into a sum of sinvsoids 
cf different frequencies. The pictorial refpresentaticn of 
tke Fourier transform is a diagram which displays the 
amplitude anc frequency of each of the determined sinuscids. 


Mathematically, the AMEFLITUDE SPECTRAL CENSITY cr more 
generally tke Fourier transform of x(t) is given by the 


relaticnshif 


-j2nft 
at) = | x(t)e dt (t=) 


where x(t) is the waveform to be decompcsed into a sum of 
sinuscids, X(f) is the Fourier transform of x(t), and j=\-1. 
The INVEESE FOURIER TRANSFORM is given by the relaticnship 


oo 
+ 2a © 
x(t) = X.(Z)e : df (1-2) 


=00 
The Fourier transform and its inverse may conveniently be 


written as 


v(t) ae V (f) 


FOURIER TRANSFCRM FAIR 
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For discrete signal processing the equaticns take cn the 
fcllewinc fcrm. For a periodic 


functicn x(t) the 
transfcra is defined by 


Fcurier 


—j(27/N) nk 
X¥(k) = X(n)je 7 K=0,, lence NS (i=3) 


n=0 


The "true" frequency component is given by 


N-1 


-j(27/N) nk 
X(kyNT) = )x(nT)e SA me (1-4) 


n=0 


efsef'-5; and f£1-&] are 


knewn as the JCISCRETE FCURIER 
TRANSFORM (CFI). 


The inverse DFT is given Ly the 


fcllowing 
relaticn 
N-1 
+j(27/N) nk 
x(n) = (1/N) X(k)je ,O=05 laces phn (1-5) 
k=0 
ard for "true" time 
N-1 
+5 (27/N) nk 
x(rT) = (1/N)) X(kK/NI)e (1-6) 


k=0 
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where 


N —> total number of samples cf tite 
and samples of frequency 


— —> sample period 

n -> sample period index 

k —> harmonic index 

F = (1/NT) —> fundamental frequency 


Fer future reference, the basic preperties cf the 
continuous Fcurier transform and discrete Feurier trarsform 
are Summarized in Table-I. Fige. gives a fictorial 
representaticn of some of the most commonly used fFcurier 
transfcrm fairs. 
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The twe gpretlems most often encountered in usirg the 
discrete Fcetrier transform are aliasing and leakage. 


1. Aliasing 


Aliasirg refers to the fact that the high-frequency 
ccmEecnents er a time function can impersonate low 
freguencies if the sarpling rate is toc low. This prceblen 
1s ccrrected by sampling the signal at a rate 2 the Nyquist 
Freguency, i.e., cemanding that the sampling rate Ee high 
encugh fcr the highest freguency present to be sampled at 
east fwice during the period. Fig.{2] gives a fictorial 
representaticn of a signal without aliasing. 


(a) h(t) (d) it |H(2)| 





cc) 














1/The— if 


(f) le( 2) #a(£)| 





Figure z - Fourier Transfcrm of a Wavefcri Sampled at the 


Nyguist Sampling Rate. 





wken the sampling period T of Fig.[{Z-c} is increased 


e= cmowm ge fig.{3-c.], the equidistant ispulses of (f) 
beccme ecre closely spaced, Fig.[3-d]. Eecause of the 
decreased spacing of the frequency impulses, their 


ccnvcluticn with the frequency function, H(f), Fig.[{3-b 
results in tke overlapping waveform illustrated results in 
the overlapping waveform illustrated in Fig.[{3-f]. This 
overlap is tke result of aliasing by the high freguency 
cc&pcnents. 


The eliasirg phenomencn will occur anytime a sampled 
analog signal contains fréeguency components greater than 
half the saapling freguency, (f )/2. The aliasing frecuency 

s 


relaticnshig is such that ali frequencies higher than tne 
analysis ktancwidth (f£ = (£ )/2 ie “tcldt = dewn jt ometha: 
max Ss 


band in an accordian-like manner with a fold at every 


Mmotipie cf f ; 





lax 
(a) n(+) fo) Alea 
7c 
(c) 4 Q(t) (d) One) 









t 
peo 


(£) a (2) *O( £)| 





Figure 2 - Fourier Transfcrm of a Wavefcrm Sampled at Less 


Than the Nyguist Sampling Rate. 


N 
ta) 








Ze Leakace 


If a periodic , band limited functicn is sampled and 
truncated tc ccnsist of other than an integer multiple of 
the pericd, the resulting discrete and ccntinuous Fcurier 
transforges till differ considerably. This proklenr is 
inherent in the Fourier analysis of any finite reccrd of 
cata. The reccrd tas been formed by looking at the actual 
Signal fer TI seconds and ky neglecting everythince that 
happened tefcre cr after this period. As shcwn in 
Fig.{4-cj, this is eguivalent to multiplying the signal by a 
rectangular tindcw. 


If the ccntinuous Fourier transfcrm of the fure 
ceSine wave cf Fig.{4-a] had been found, its contrikution 
would have teen limited to a pair of impulse functicns cn 


the frequency axis, Big.{4-a}. 
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Ficure 4 - DFT of a Eand-limited Feriodic 


Wavefcrm:Truncation Equal te an Integer Multiple. 


Zc 








Hee eecrfiect cf truncaticn at other than a muiticle of 
the fpericd cf tke sampled signal is to create a fericdic 
function swith sharf discontinuities as shcwn in Fig.[5]. 
Tre muiltiplicaticn by the data window in the time dcmain is 
eguivalent tc performing ccnvclution in the frequency 
demain. Ccnseguently, the frequency functicn is ne lcnger a 
single impulse but rather a continuous function of frecuency 
with a lccal saxinuz centered at the criginal impulse anda 
series cf sftrious peaks termed sidelobes. These sidelches 
are respensitle for the additicnal frequency ccompecnents 
Nleakage", which cccur after frequency domain sampling. An 
expanded view cf sampling in the freguency domain is shown 
eee ig.{ G}. 
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peogure © - CFI cf a Periodic Waveform:Truncation Interval 


Not Fgual to an Integer Multifle 
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H(f)+Do(t)ety SeeTe) 





Figure 6 - Expanded Illustration cf the Convoluticn of 


Figure [{5e]. 


Tre tsual appreach tc reduce the agrcunt of "leakage" 
thrcugh tkese sidelobes consists cf applying a data window 
cr weighting functicn to the time series, which has lcwer 
Sidelobes in the freguency domain than the rectangular data 
window. Tle process Of weighting consists On 
agplitude-scculatirg the signal prior to Feurier analysis. 


[ie welighbtinc ¢recess is illustrated in Fig.[7]. 







FOURIER TRANSFORM 


eye 
TIME OF PRODUCT EQUALS 


FUNCTION FOURIER 


TRANSFORM 





TRANSFORM OF THE 
TIME FUNCTION 
CONVOLVED WITH 


WEIGHTING TRANSFORM OF THE 
FUNCTION WEIGHTING FUNCTION 
Ficure 7 - ‘The Weighting Functicn Frocess. 





Ikere are many such windows but one of the SIPCLLSE 
and most cften used is the Hanning function. This function 
feed CCSine-rell curve given by 


w(t) = 1 - Iecs2qt Q<t<TIc — 
ae eae Aas!) 


where IT 1s the truncation interval. The magnitude of the 
ce 


Feurier transform cf the Hanning functicn is given by 
{mrt} | = I0(f) + 1 £+1 + an 8 
= (z) ql o¢ qe! C to) 3 ( ) 


meere C(f) = (Sin("Tef) jyr£t 


The Hantirg function Fourier transform fair is shewn in 
pac .{ £j. 


Windcwirg of the input frame by the Hanning furcticn 
recuces "leakage" Ey imposing a gquaSi-pericdicity cr the 
input sicral and eliminating the discontinuities Letween 
periods. Tke Hanning function widens the main lobe of the 
(Sin x),yx shape and reduces the spectral amplitude Ly a 
little acre than 4dEB, but it causes the “interfering" 
e¢idelcbes to fall off at 18d3 fer octave instead cf €dE as 


mrethe crectarcular window. 
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AMPLITUDE 


Change of Scale 





.02 


(b) 


Figure € - Eanning Functicn Fourier Transform Pair. 








It shcula cre remembered that the non-zero frequency 
ccaponents are considerably EROADENED cr SMEARED with 
respect tc tke igpulse function. In general, the mcre the 
leakage is reduced, the broader or more smeared the results 
of the discrete Fcurier transform appear. MTable-II cortains 
ecme cf the acre ccmmcnly used data windows. 
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TABLE-IT 
Data Windows 


3-dB NOISE HIGHEST ASYMPTOTIC 
WEIGHTING FUNCTION BW Bw | 2LDELOBE ROLL-CFF 
(dB/OCTAVE) 


RECTANGULAR 

0.85B/ 1B 
TRIANGLE 

e255 iL. 358 
COSINE 

1.273 |P, 205 


HANNING = (COSINE)@ 
1.48] 1.5B 
(COSLUuE) 3 
1.61B |1.73B 
(COSINE) 
1.888] 1.9B 
HAMMING = (COSINE)© + 
3% PEDESTAL 643/0CT 
1.3B/1.36B 3EYOND 
5B 


ge 








DPD. SFECTRUE ANALYZER PERFORMANCE FIGURES OF MERIT 


1. SasgEling Frequency and Anti<Aliasing Filter 


Ike sauplirg freguency must satisfy the Nyquist 
Sagpling tkecrem cr aliasing will result. Fer data that has 
conpcnents Eeycnd the analysis range, the componerts are 
attenuated Fry means of a lowpass filter, known as an 
Panti-~aliasirg" filter. Awe typical “fall Wort" irat< sror 
anti-aiiasine filters in today's equipment is 120 dByoctave. 
In practice, a Sampling rate of approximately three times 
the width cf the desired analysis range is a gocd 
CcBEromise. Sampling at higher rates requires more 
equiprent, thich is undersiratle fcr eccncmic Lreascns. 
Sampling at lcwer rates fails. to take “into acccunt the 
actual fall-cff of available anti-aliasing filters, anc thus 


dees not e€lininate aliasing ccnpletely. 


w~ 


2. Eaadwidth (Freguency Range) 


{ke kandwicth determining parameters of a spectrum 


analyzer are: 


1. The memory capacity (which determines the 
curker of Satples availaktle fcr prccessing). 


Brling frequency (which, in combination 
witht 1, détermines the length of signal (ir 
s) which is available fcr prcecessing). 


3. The shape of the "windcw function" also kucwn 
as the "weighting functicn". 


dad 
la 





3. &escluticn ce 3-dB Bandwidth 


om we a ee ee eee 4 — oe = aw ee ae eee 


in Takle-II the 3-d5 bandwidth for various weighting 
furcticns were ccmpared using the relation 


B(H2) = 1 


a a SS CESSES 


heighting-functicn length (sec) 


The 3-dE FEandwidth is a good measure cf a spectrum 
analyzer's akility to resolve two equal-amplitude ad-acent 
sire waves. Studying Table-II, it apfears that the 
rectangular weighting function has the best resoluticr, but 
Momeceality, the respcnse characteristic is limited in 
Meefvulncss Ey the side-lobe structure. In practice, it is 
mcre impcrtant to resolve unequal than equal anplitude 
adjacent fréeguency components. 


4 Neizse Eandwidth 


Scmétines called “effective bandwidth", it is used 
to ccnvert cr normalize power spectrum measurements to power 
spectral censity (gcwer per unit frequency). The ncise 
bandwidth is tke bandwidtn cf a hypcthetical rectangular 
filter, which passes a signal with the same mean-square 
value as the actual filter when the filter input is white 


Gaussian rcoise. 


Se eebiguest Sidclobe Wevel an@ Asymptotic Rcli-ott 





It is useful to specify these parameter, since it 
prcvides scme indication of the prceminence cf the side-lohbe 


structure. The largest sidelore is generally the first cne, 
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with the Hamaing functicn being amcng the exceptions. The 
third side lcte of the Hamming function is the largest. 


Ike presence cf side lotes makes it difficult te 
distinguish -Letween a high-level frequency component and a 
lew-level ccufenent which is close in frequency. Fer this 
reascn, lLectangular weighting, which leads to a spectrum 
response eégual tc the true Fourier transform of the 
tigre-truncated Signal, is used only cccasionally when 
analyzing trarsients cr shocks. 


angis-Sided Spectra 


ON 
e 
tn 


AltrEcugh it is nct an actual measure cf performace, 
a discussion of single-sided spectra is warranted at this 
time so thet the results cbtained fren the Varicus 
algcrithrs is fully understood. 


Ike firite discrete transform presumes that the 
ingut wavefcrm is the sum cf ccsine-phase and sineé-phase 
frequency ccuafenents, all harmecnically-related to the lowest 
freguency tkat can be recognized. That fundamental is the 
freguency thet cempletes one cycle in the length of the 
input frame. The calculated harmonics range fren ‘d-c) 
through 1 (tke fundamental),2,3,etc., to Nyz, the tGaximnun 
reccgnizaktle karmcnic, which completes one cycle in tke two 
mmecE Val: associated with twe pents of the input (1.€., two 


Sempples per cycle). 


Ir tke general case, the discrete Fcurier transform 
Can cperate cn a complex time series input functicn and 
prceduce a ccnoplex spectrum. In such a case, an input 
furcticn cf N ccmplex values (i.e., 2N independent real and 
igeginary farts) yields a spectrum of 2N independent real 


and isegirary farts cepresenting N cc#rplex frequency 
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ccefficients centered about 0(d-c) and ranging from harucnic 
=—872 to +N yz. The "negative-frequency" ccmponents are an 
integral and necessary part cf the function. 


A real input function, as an electrical waveforn, 
may te ccrsicered fcr this general case as a ccmplex 
furcticn whcse isaginary parts are all zere. The discrete 
transfcra cf such an input yisclds a complex spectrum as 
accve. For this special-case input however, the spectrum 
will ke syzametrical akcut the center or d-c point. The 
necative-frecuency parts from 0 to -Ny2 are an image cf the 
Ecsitive frequencies, the real or coSine parts having even 
symoetry ard the lmaginary, sine Farts odd. The 
necative-frecuency half thus becomes redundant, and an input 
of N real feints yields N independent reai and inacinary 
parts cf Ny2 ccemplex spectral coefficients. 


Ncnetkeless, the negative-frequency coefficients are 
still a necessary part of the spectrur in that they 
represent kalf the amplitude of the harmcnic components. In 
cther words, if the input functicn contains a cosine-phase 


ccgrpcnent cf amplitude A at frequency eat ther the 
c 

resultant pesitive-frequency ccefficient for frequency 3f 

will have agplitude A/2 and the negative-frecuency 


th 
@eetticiert fcr —, will also have amplitude Ay2. The 0 


or d-c ccefficiert, however, dividing cr being “skared" 


ketween tke pesitive and negative halves cf the spectrum, 
Hii have ft1l amplitude A for an input d-c level of A, 


Meaeh is tte gciat cf this discussion. 
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Pee sPRCIRAL ANALYSIS ALGORITHMS 


Ite fast Fourier transform (FFT) is a Eighly 
efficient frecedure fcr computing the DFT cf a time series. 
It takes advantage cf the fact that the calculation of the 
Seetcticients cf the DFE can he carried cut iteratively, 
which results in a considerable savings of ccmputaticn time. 
Tc gain insicht to the FFT algorithm, a general develorfement 
is freserted fcllcwed by several examples. A more in depth 
Seudy of the FFI may be found in Refs.[6 i, [16], and [{ 17]. 


a. General Cevelopement cf the FFT Algoritho 


{ke discrete Fourier transform Eq.[1-3]j is 


N-1 
=j (27 /N) nk 
X(k) = x (n)e pok= Ogu liarctecety Wee (=3)) 


n=0 


which gay ke written as 


N-1 
nk 
X(k) = x(n) W (1-9) 


n=0 
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wEere 


=n 7 N 
Ww =e (1-10) 
| pk | 
As can be seen, W 1S Peblodie sand) Of e per iced en. 
eves , 
(n+mN) (k+1N) _nk 
” = W ev m,l= OF seis (111) 


In carrying cut the FFT algorithm, the symmetry and 
a « e nk s e 
Perdodicity cf i are exploited to achieve an increase in 


efficiency. Tc emphasize this periodicity, W will be 


resrittec as an in carrying cut the general developemert. 


nk 
od.) 1-3), w is a complex exponential, thus 


an €xaBiraticr of £Fq.,1-9] shows that, in the case when x(n) 
is a ccmgplex sequence, a complete direct evaluation cf an 
N-point CFI reguires (N-1)2 complex multiplicaticrs and 
Mam—1) ceBplez additions.! fFig.{9] shews the compariscn of 
multiplicaticns reguired by the direct calculaticn cf the 
mete and fcr the kase 2 FFT algcrithn. 


1In mest literature, N is assumed to be large, 
sc that (N-1) can ke approximated by N. 
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Direct Calculation 
{ 







256 


NUMBER OF MULTIPLICATIONS (x 1000) 






128 


Al 


64128 256 512 1024 
N {number of sample points) 





Figure 9 - Multiplications Required to Ccmpute the Dirt. 
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Es Eecimation-In-Time (DIT) Algorithn 


The general principle tehind the FFT is to break 
tke ocriginal N-point sequence into twce sherter sequences, 
the DFfT's cf which can be combined to give the DFT of the 
criginal N-pcint sequence. 


Assuming N iS a power of iZ, define two 
(N72) -poirnt seqtences en) and hue as the even and odd 


meakers cf tke seguence x(n) respectively, i.e., 


x, (n) = x(2n) n= eda - 1 
Cig=12 ) 
fn) = x(2n+1) n= Dot o an - 1 
2 2 
The N-point DFT of the sequence x(n) can be written as 
N-1 N-1 
X(k) = x(n) i - ) x(n) (1-13) 
n=0 n=0 
n even eOec 
N/2-1 N/2-1 
= x(2n) We + x(2nb1) Wyott) (1-14) 
n= n= | 
Then using the following identity 
o Bea ie 207 (N72) : 
Eq. [1-14] Can written in the form 
N/2-1 N/2-1 
yok nk i 
x(k) = x4 (n) Wy 7, + X_(n) Wy 7. (1-16) 
n= no 


4Q 





SSAC) + Wako (K) (127) 


Each of the sugs in Eq.{[1-16} is an Ny2-peint 
me, the first sum being the N/2 point DFT cé£ the 
even-nuagterec epeirts of x(n) and the seccnd peing tke N/2 
pemee, CEI cf the cdd-numbered points cf the criginal 
seguence. Zlthcugh the index k ranges over N values, 
K=C,1,...,N-1, each of the sums need only te computed for k 
Eetween 0 arc N/z-1, since at) and oD are each fpericdic 


mm With pericd Ny2. Aiter the two DFT's ccrrespondirg to 
the two sus in £Eq.{1-16] are computed, they are then 
ccgtined tc yield the N-point DFT, X(k) as given by 
moar 1-17). 
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2. Ike Chite-2-Irarsform (CZT) 


TEe chirpe-z-transform (CZT) is an algoritkm for 
evaluatin¢c tke z-transform of a finite duration seguence 
alcng certain general contours in the 2z-plane. The 
Ecilcwing derivation is based upon a regpert by Rabiner, 
schater,ard kader [18]. 


ae Ike Cerivation of the CZT Algorithno 


In gereral, the z-transfcrm cf a seguerce of 
Samples {x(n)} is defined as 


X(2) = x01) 52 7 (1-18) 


Assummin¢g tkat the right side of Eg.{1-18} converges for 
some values cf z, and restricting the evaluation of the 


z-transfcro to a finite duraticn sequence of 
Sairples, pea d= 10 |. 
Ne | 
al 
X(z) = “%(n) Z (1-19) 
n=0 


The special case cf the z-transform which has 
received ccrsideratle attention is the set cf points egually 


Spaced arcund the unit circle of the z-flane, i.e., 
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a = exp (j2TK/N), k=0,,1j0%0,N=1 (1-20) 
which gives 


N-1 
-j}(27/N) nk 
X¥(z) = x(n)e p K=O pecs N (1-21) 


n=0 


By ccmparing §Fg.{1-21] and Fq.{1-3], cne can see that it is 
the discrete Feurier transform (DFT). The DFT can be 
evaluated scre efficiently by using the FFT rather thar the 
CZ1T, but tte latter allows cne to evaluate Eq.!-19]j along 
the mcre cereral ccntcur 


Z = Ad Pie <0 caer) aa (=22) 


where M is ar arkitrary integer (not necessarily equal tc N) 
and A and W are arkitrary complex numbers defined Ey 


to» 
ul 


A exp (i2Te ) (1-23) 
o ° 


= 
| 


HW exp (52nd ) (1-24) 
° ° 


Fig.{10] shcws the significance of A , et 2 and 2. in 
OQ 


the z~plane. 





Z-PLANE 
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gre 10 > Contour of the CZT in the 2-Planc- 
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me: the Specialcase of the DEI, A=1,) B= 3, and 
W=exp(-jz7m/sN); thus the DFT coefficients of a finite 
duration secuence are the values of the z-transform of that 
same segvuence at N evenly spaced points around the unit 
curcle. 


The parameter W determines the rate at which 
Oo 


mien CONtGCYr spirals: if W& is greater than unity, the 
O 


ccntour spirals toward the origin as K increases, and if id 
Oo 


1s less thér unity, the contcur spirals outward as K 

increases. TIte parameters A and 2790 are the locaticn in 
O o 

radius and angle, respectively, cf the first sample, i.e., 


fer k=0. Tke remaining samples are located aicng the spiral 
Centcur «itt an angular spacing of 27%. The equivalent 
Cc 


e-riam@e ccntcur of Frig.{[10] is shown in Fig.{11]. 





a 
on = min Ay 


Figure 11 - Contour of the CZT in the S-Plane. 





Tké equivaient s-plane contour teginning foint 
may be fcund ty making the substitution 


syd 
Z=e 
tken 
2 = A = exp(s T) 
G oO 
cr 
Sem =) (171) 1lpAa 
e 
=eoemt JW = (1/7) {inA + jJ2778 ) i= 25) 
Cc oO oO oO 


Pare d geteral] point cn the s-plane contour 


1 


exr(s TIT) = AW 
E ( k ) 


N 
it 


tren 


fn 
" 


=K 
Cer 2a (AN. ) 


anc finally 


in 
i 


(1/1) (Lna-klowW) = s -(kyvT)low 
Cc 


Mrge the Eefiraling contours in the z-plane correspend to 
straight lines in the s-plane. The rate cf spiraling in the 
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z-Flane cetermines the slope cf the lire in the s-fplane. 


The task now at hand is to ccmrute 


Eq 12] 


aleng a general ccntour given by Eq.{ 1-221. The evaluation 
of 2(z) at Pa will be written using the notation xX, 
k 


Le 


cr 


where N is tke length of the sequence x(n). 


Making the follewing substitution, 


Sradgnanated ky Eluestein [5 }. 


nk = n2 ¢ k2 - (k-n)2 
2 


gives 
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(1-27) 


(1-28) 


whick was 


(ii=2 9) 





N-1 
See nae 


<< 
Nt 


x(n) A 


n=0 


Ne 


SNe cash oy 2. — k= ne io 
Peay cwe) iw ae ee 


(1-30) 
n=0 

Althcugh fFq.[{1-30] appears to be mere ccmplicateé than 

Eg.{ 1-28], it allows efficient hardware isplementaticn. 


Eq.{1-30} ray te viewed as ktasically a three step precess: 


1. ire flEsSt step ccnsists (ie a ccuplex 
pre-ultirflication forming a new sequence 


“Dn ne@/z 
¥te) = x(npA WwW pe 0 lg sey NO (1-31) 


2. Ire seguence y{a) is then ccnvclved witk the 
segqtence v(n) defined as 


~n2@/2 
v(r) = i fi (1-32) 


tc give the sequence g({k) defined as 


N-1 
g{k) = y (n) v (k-n) , K=O 4 Aig sey te (1-33) 
n=0 
ae k2/2 
3. The sequence g(k) is then post-multiplied ty W 


tc cive 
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N-1 


k2/2 
X =wW y(n) v (k-n) ~k=0,1,...;,8-1 (1-34) 


n=0 


A pictcrial representation of the processing operaticns of 
imemcal algcrithe is shcwn in Fig.{12]}. 


n CONVOLUTION 


zZ 
* woh UP 
vee 





Z Z 
AT Dyn ie we je 


Figure 1Z - Precessing Operations of the CZT Algcrithn. 
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The adventages cf the CZT algorithm cver the standard 
FFI are Sugnarized below: 


1. oN, the number of points in the input sequence, 
dcesr't have to be equal to M, the number ot 


4 
? 
points é¢t which the z-transform is evaluated. 


2. Weitber N nor M has te be composite - beth can 
re -rime numbers. 


aie The Sear eang point in the z-plane and the 
angular ¢facirg of the ie are arbitrary. Thus 
the frequency resolution and the range of 
fregtencies are arbitrary. 

Hs dre contour does not have to be a circle iz 


the z-flane. The CZT algecrithm can evaluate the 
z-transfcro along spiral contours inside anc 


Olteide the wnit circle of the z-flane. By 
ean wating the z-~transforn oe a selected 
ceontctr, the transfer functicn of a linear syster 


can re ad=usted to sharpen or_ flatten the 
£regquency resfonse. The cheice of centour is 
depéndert on which of the system's poles should be 
emphasized. 


, SSN, W= Sea ee then the CZT can 
evaluate the DFT, even when N 1S prime. 


6. Rith the CZT algorithm, interpolaticn between 
time samples cf a andlinited functicn can te 
CkEtained with greater ease with the DFT, and 
interrclaticn at arbitrary sampling intervals Can 
be perfcrmed efficiently. 
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@. ‘IFe Prime Transform Derivation 


The foll¢gwing derivaticn is based on a paper by 
Whitehouse, Means, and Speiser [26]. A brief review cf the 
nugkter tkeory properties used in this derivation is given in 
Appendix B. | 


In the specific case when N is an odd frime, 
Bq.{1-3] can ke written as 


N-1 


Pic) = ot (1=35) 


n=90 


tc give tke (D-C) frequency compcnent and 


N-1 
sean k /N 
X(k)-x (0) = x(n)eé pKa lige ule gat (1-36) 


a=1 
te give tke remaining components. 

Since cnly non-zero values of n and k cccur in 
tke right band side of £q.{1-36], it is possikle tc replace 
the product rk by @ primitive root raised tc a sum [esi 

Preceeding , in the summation cf Eq.{ 1-36], the 


terms are permuted and the order of the eguaticn changed by 


om 





the fcllcwine transfcrmations: 


(n)p —= ((r69)P)) = r’ modulo N 
- (37) 
r- modulo N 

gly <b ae auiicoal 


where "©" is a primitive root and "pf" isplies permuted. 


()p —— ((r'*)P)) 


N-1 0 
Mcting that rc modulo N = cf mcdulo N, Egq.{1-36]}] is 


written as 


N=-1 
: = ,2m_.((n) p+(K)p) 
F (C802) =o =) HD app, RSID) 
Se (1-38) 


An examination of Eq.| 1-38], shows that the sequence 


cr 


we ()p) 7 x(o) Pomthiewncircular Cormelatlonemon the 
r 


sequence (exp -aten/ny eo) ‘ 


In @atriz ncetation, Eq.[ 1-36] may ke written as 


X - x(0)D = F x (n) (1-39) 


wrKere 


D= |. (1-40) 





' ' 
and X ard x (n) are cclumn vectcrs of size (N-1) derived 


frem X(k) ance x(n) by deleting, respectively, X(0) and x(0). 


q 
The wWatrix F Can be factored into three 


matcices 


t 
fe PCE (1-41) 


where E is an (N=1)  permutaticn matrix, € 94s ap 
t 
{ (N-1)X(N-1) ) circulant matrix, and P is the transpose of 


iP. 


The elements of the C matrix are 


ce = F De k=1. eal 1-42 
r yk (anieeyp 


where "r" is a primitive root cf N and "p" ipplies permuted. 
Tie ele@exzts cf the permutaticn matrix,P, are 


= § Aipk= Viger, Non (1-43) 
n,k (n)p,k 


where the <«kronecher delta functicn. The 


40) Eek 
circulant matrix,C, can be igplemented with a recirculating 


trérsversal filter whcse tap weights are determined Ey 


h(x) = y nF , N=1,...,N-1 (1-44) 
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Thus the prime transform algcrithm offers more 
sisplificaticr than the CZT algorithm because in hardware 
isplementaticn it is possible to eliminate the pre- and 
pest-multipliers. The processing operations may be 
sum@arized as 

1. pergottaticn cf the input data 


peeeeitretriatr correlation 


3. permautaticn of the output Fourier ccefficients 
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F. HARD RZRE AND SCFTWARE IMELEMENTATICN 


in kEcth tkardware and software implementation of theses 
algcrithas, the gcal is to have the ability tn do real-time 
Signal frecessing fcr the frequency range cf interest. In 
spectral anélysis applications, signal frosessing takes 
fFlace in real] time when the spectrum level cf each frequency 
ccmponent is updated at a rate equal te at least the 
tandwidth 6 (82) cr the systen. 


Cependince upon the system, the signal for processing may 
Me in any cf tktree forms. 1.¢€., 


1. Analog-> toth time and signal amplitude are 
CclLtisanuous 


2. oes analog-> time is discrete and signal 
agplditude is continucus. 


3. Dicital-> time and signal amplitude are tctk 
discrete 


oan 
fos 
ry 
o 
f= 
a 
a 
Iw 
ip 
a 
\+ 
(4 
I 
tt 
(e 
D 
im 
tr 
be: 
ttn 
{po 
ler 


=~ ee ee oe 


a. Scanning Analyzers 


Ike conventional scanning-filter analyzer 
ccnsists cf a single bandpass filter, whose center frequency 
is “shifted" by using freguency-translaticn techniques. 
Typically a single bandpass fiiter is employed and the 


filter analyzes cne frequency at a time. 


= 





Since a single filter is employed and is afrlied 
successively tc each frequency, it must dwell for a period 
of i/B8 at €ach frequency. Thus (14)X Ne. of sters equals 
tke tctal analyzing time "T", The narrower the filter the 
lenger is the tctal sweep time. Since each frequency 
ccipenent is upcated only once per time T, real time signal 
precessirg is ret fossible. 


FE. kultifilter Spectrum Analyzers 


Ccrstant Fandwidth multifilter analyzers were 
develcped fcr real time analysis. In this tyre of analyzer, 
the Sicnal passes through many parallel filters 
Simtultanectsly , typically 500 contiguous magnetostrictive 
filters fer analysis range. The filters and the asscciated 
circuits gust fossess uniform and constant gain so that a 
fixed amplitude sine wave of any frequency in the analysis 
range results in a spectrum output which is ccnstant. 
hcwever, frecuency flatness of better than +3dB is difficult 
tc maintain, thus the magnetcstrictive multifilter analyzers 


are net used extensively. 


2. Ligital spectral Analysis 


= 


tke 611 digital spectrum analyzers employ a computer 
cr ccgputer-Jike circuitry which is programmed to inplement 
Meee brit alccrithm. Some of the structural factors fcr the 


desigr cf such a system are discussed below. 
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a. Technolcgy 


Easic speed and functional biccks available are 
dictated Ly technclogy. Table-III shows a ccmuparison of the 
ajcr legic fagsilies. Because Of speed-ccmplexity 
tradeoffs, tte cheice of logic family is usually not a 
clear-cut decesion. This complexity encompasses many aréas, 
but important are available MSI and LSI functional ucdules, 
system interccnnections, multiple function mcdules, and type 
of multiplier. Tc gain speed in the FFT algcrithm a fast 
Bietiplier is reguired in the butterflies, thus the type of 
multiplier is an inportant consideration. Table-IV is a 


ccgsfariscno of zultipliers. 


TABLE-IV 


MULTIPLIER COMPARISON TABLE 


Computation Time | 
Type of | Component (nsec) | Package Count 


Configuration| Basic Package | Technology |16 x 12 16x 8/9 x91 16x 12V 16 x 819x9 

















Fast trapezoid | Custom two- ECL 32 28 | 28 He Min) 82 1) 45 
(Fig. 8.35) bit adder 
package 
Tree Commercial EGIL 62 48 40 160 120 85 
four-bit 
adder | 
package | 
Tree Commercial TLE a 100 85 160 120 85 
four-bit : | 
adder 
package 
(2 x 4)-bit Tre 261 195 135 80 58 30 
(commercial) 
array pack- 
age 
Two bits at a Custom two- EGL 130 120 100 40 36 25 
time add- bit adder | 
shift package 
Two bitsata | Commercial | EGE 210 | 140 116 40 SG) \Aieos 
time add- two-bit 
shift adder | 
package - | | 
Two bits ata Commercial | ThE 450 300 250 | 40 36 25 
time add- two-bit | | 
shift adder | | | 
package | 


Taken from Ref.{17] . 
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Ee 2igoritha 


Ihe structure of the FFT algorithm is a major 
factor in detergining whether or not special-purpose 
hardware cr a general-purpose computer is used in a 
particular apelication. Table-V¥ lists the aritkmetic 
oferaticnos recuired for different radix algorithms. When 
isplemented with software, the radix-8 algorithm is near 
Cptimum. dHcwever, it lacks the flexability required of a 
General-furpese ccnputer and thus is only used in sfecial 
PUtEcCse arpriicaticns. Although the radix-2 reguires 
additicnal ccrputations, its simplicity has made it foflular 
for use iz general furfose computers were real time analysis 


is ret a rectirement. 
c. Eardware Architecture 


then insplemented with hardware, the radix-4 
algcritha hes received the anost attention. Hardware 
igplementaticr cf the PPT algcrithm has gained speed by 
trading cff flexability. The speed is gained by a ccmplex 
nix cf severai types of parallisn. One such type is 
pipelininc, i.e., the process is broken uf into several 
sequential tasks and executicn proceeds assenblv-line style. 
Ey pipelining many processing steps can be progress at any 
given tine end thus perfcrming real time analysis. Given 
meat an B-pcoint FFI requires eo butterflies, there 


are four kasic ways tc organize the prcecesscr. 


5g 
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1. Seécuential-> all a ogee butterflies are 
Ferfcrmed sequentially. 


2. Cascade-> if log N arithmetic units (AU) are 
° < ® 
available, the algotihm can be firelined sc 


that each AU computes (N/2 kutterflies 
befcre advancing data through the fire. 


3. Fatallel-Iterative-> if (N/2) AU's. are 
aveilable, the algorithm can ke precess in a 
log N stage pipéline. This requires oe 


tutterfly times. 


4. ALrray-> if ener AU's are available, the 


et array can pipelined into cne executicr 
ime. 


Tarle-VI ccufares these different organizations fcr an 
N=1024 PFT arc an execution time of 1s The total execution 


tigze does nct take into account the I/O overhead tine. 


TABLE-VI 


Comparison of Processor Organizations 











Machine . peach. Unite Rota 
Organization (Butterflies) Execution 
Time (ps) 

sequential tL Buz, 

Cascade LO, Sie. 

Parallel 
Iterative pics 10 
meray 5120. de: 
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In addticn tc fipelining, among the tyres cf parallism are 
1. Tine cverlap cf control and memory functicns. 


2. The addition of a small amount of high speed 
LENCry. 


3. Eicher radix algorithms. 


TaEle-VII ccapares some isgportant farameters of a few 
modern sicgral] precessing systems. 


Like the FFT algorithm the CZT algcrithm can ccmrute 
the frequercy ccmpcnents of a signal mcte efficiently than 
by direct evaluaticn of the DFT. However, for relatively 
large values of N, the CZT algorithn is less efficient when 
Gc@raread tc the FFT algorithm. For this reason interest to 
digitally isflement the CZT algorithm waned while great 
strides sere being taken to develcpe FFI software and 


bharcware. 


The GZT does offer many Signal prccessing 
advantages, Evt cnly if it can be implemented using sampled 
analog teckriques. With the developement of Surface 
Accustic wave (SAk) and Charge fransfer Devices (CTD), 
interest in the cCZT has been re-kindled. The prime 
advadvantace cf these devices cver digtal equipment is that 
when the ccrtinucus analog signal is sampled the analog 
ides sis retained, while in a digital system the analog 


value must Le guantized which is a source cf error. 


CID¢s, which include the Charge Coupled Device (CCD) 


ard Eucket-Erigade Device (8BD), sample the Gpcical OF 


OZ 





electrical infut and convert it to a packet cf charge whose 
magnitude is prcfortinal tc the input. The charge racket 
moves down the CTD stages, shift register fashion, tunder the 
Ceatrcl cf clock voltages, The charge packets are 
transferred fren one stage to the next until they readh an 
cutput stace where the magnitude of the charge facket is 
sensed anc ccnverted to an electrical anlaog signal. The 
cleck freqtency and the number of stages transferred 
determine tke time delay introduced. 


In é@ SAW dévice a mechanical ripple, cr wave, 
travels alone a sclid polished surface in much the same 
Manner aS an ccéan wave travels through water. Altzinun 
electredes are usually deposited on a fiezcelectric material 
such as STI guartz, and they act to launch and receive the 
accustic saves. These transducers convert energy bketween 


electrical ard mecharical domains. 


In a CCD the delay line is tapped and weighted by 
using a split-gate weighting technigue which is inccrtecrated 
intc the mask at the time of fabricaticn or ty weighting the 
taps witk external resistors. The tap weight ina SAW 
device is photc-lithographically determined and is 
MEerperticral to the length of the interdigitaticn which 
accepts tke accustic wave. The pre-and post-multigflies are 
acccmplisked usinc wideband four analog aultipliers with 
weights generated ky programmable read orly wemories (FCM) 
and ccupled thrcugh D/A's or by multiplying D/A's with FROM 


weight stcrace. 


Ike transversal filter implementaticn of the C2T 
alccritha achieves its high speed not by reducing the rumber 
cf multiplies like the FFT, but by putting the equivalent of 
M earalle] evltipliers all on the same chif. This reduces 
tthe number cf ccmputational steps fron odes for the FFI tc 
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merce tue CZ1T. Alsc the N additons per component are 
Sinultanecus. Currently, the CCD CZT can operate in real 
tige at speeds up to SMHz and the EBD CZT is limited tc a 
sample rate cf a few hundred KHz. For higher sample rates, 
up tc several FHZ, SAW devices are attractive. [In addition, 
tke sagfled analcg devices require low pewer and are light 
weight. The advantages of the digital FFT are its 
fFlexability and accuracy. Programability is kteing develcped 
meee the €CD and SAW filters with up te 128 digitally 
Switched taps have been built. 


As rected earlier, the prime transform cffers 
aoditicna ] Simplicity by elirininating the multifliers 
reguired in the CZI implementation. The multipliers are 


replaced with analog permutéer memories. The commercially 
avéilable sé€rial access memcry stores analcg samples as 
charges in an array of MOS capacitors under the ccntrcl of 
read in ard read out shift registers. Currently the maxinun 


sarple rate for this type of device is 5 MHz. 
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i 
(LIT) ALGCRIIEM 


ZO abluctrate 


peints of x(n) is assumed to Le a fower 


the 


In this example N=4=2e and ¥ =2. 


For N=4, Eq.{1-3] can be written 


X(0) = x(0)w° + 
X(1) = x(o)w°? + 
X(2) = x(o)w? + 
X(3) = x(o)we + 


Se written in matrix 


X(0) 
X(1) 
G2) 
X(3) 


x(1) wo 
x(1)w> 
x(1)w* 
x(1)W? 


form 


+ x(2)w® 
+ x(2)w 
+ x(2)w? 


+ x(2)we 


= 


= = = = 
Oo £ NY O 
= = = 

<D” Ox, AD 


66 


MATHER2TICAL ILLUSTRATION OF THE DECIMATICN-IN-TIME 


FFT algorithm, the number cf sample 


¥ 


of 23 i.e. N = Z e 


+ x(3) We 
+ x(3)W? 
+ x(3)we (2-1) 


+ x(3)W? 


x(0) 
ze (222) 
cle) 


iS) 





The first step in developing the FFT algGCEiLthmr ter this 
example is tc rewrite Eq-{[2-2] using the fcllowing relation 
fics numter thecry[{23}. 


rk rk mcd {N) 
W = W (z-3) 


where {nk mucc(N) ] is the remainder upon division of nk by N. 
Hence if A\=4, 1=3,k=3 


; ==8 (2-4) 
sirce 


nk c ; y 
mee = expel (-j27/4) (9) ) = exp (-59r72) 


= €xp(-jmy<) = exp( (-j21/4) (1) ] 


1 rk macd(N 
=Ww = i oe (2-5) 


mas the cew fatrix is written as 


X(0) ior ork x(0) 
X(1) 1 ow owe wd | x(2) (2-6) 
x(2)| ta w2 ws wt | x2) 
X(3) 1 we we whl | x(3) 


= 


Mmenwseccrd sté€pf is the factcrization cf the square 
mies in Fq.f2-€]. As an aid in showing hcw Eq.{Z-€] is 
EPaetered, a réew nctation is introduced at this tine, where 


tre intecers n ané k are now represented as binary numbers. 
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a 
' 


tt 


= Gree 2g 3 > eae CU, 01, 10, 11 


K = Glens -> ar, = CO, 01, 10, e 


Tren writing the tase 10 equivalents 


De 2zdeet e kK = 2k +k (2-7) 
1 0 1 0 


Peep i=9) can new be written as 
: : (2 mache) 4 2k) 
idea + 
X(k,,k,) = a Sot ae 10 10 (2-8) 
0 oe 10 
= n,=0 


m=O 


0 


Tre aimpcrtart fpcint here, is that the single summaticn of 
Piet 9 )mecst te replaced by " " summaticns in crder to 


entrerat¢e all the kits of the binary representation of hn. 


Cepeimuirg 


Be tg) ee) *,) : eee teaou ac Po goo 


Ln k Zee (2k hn 
hi 1 ale Oa 10! 0 


2k.n (2iee eee saa 
The term in Exackets is equal to unity since 
Ln. k 4\nk P72 i lasik n.k 
ee E aes E Lol, 2der eet SS (2-10) 
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Eq.{2-8] car now te written as 


i i 
ek in (2k.+k.)n 
= Orr O20 
X(k, 5K) = ) 2 19) W (2215) 
n,=0 n,=0 
Each of the summations are treated individually anc lead 
G@itectry tc the correct factorizaticn cf Eder . 


Rewriting the bracketed summaticn 


1b 


x (KysNy) = x(n, 


na} 


1 


ny) W 


ekyny 


C212} 


Directly evaluating this summation gives 


x, (0,0) = 
x, (0,1) = 
x, (1,0) = 
x, (1,1) = 


x(0,0) + 
Berar 
x(0,0) + 


m(O,1) + 


x(1,0)w° 
x(1,1)w° 
x(1,0)W° ee 


x(1,1)W° 


re-writing in matrix notation 


© 


x, (0,0) 1 


te 


x, (0,1) 0 


S 


x, (1,0) A 


f- 


x, (1,1) | 0 


oO 
= oO 


= 
eS 


x(070) 
O51!) 

(2-14) 
SOO) 


Sek lesley) 


69 





Now writing the outer summation as 


(2k, +k))n, 


Directly evaluating the summation gives 


x (0,0) 
x, (0,1) 
X5(1,0) 


X5(1,1) 


x, (0,0) 
x, (0,0) 
x, (1,0) 
x, (1,0) 


iien in matrix notation 


x, (0,0) il: 
x, (0,1) als 
x, (1,0) 0 


x5(1,1) 0 


ONS = 
yt 


oO 
be 


- x, (0,1) W° 

2 
+ x, (0,1)W 
A x, (1,1) W> 


s ipa ie 


0 x, (0,0) 


2 


x, (0,1) 


= 


x, (1,0) 


w sen ase) 


G@e@ranine Fg<.{2-12]} and { 2-15] results in 


X(k_,k = x (k k ) 
ae 9? os a 


Ge medtircc ir @atrix form 


X(0) ae 
X(2) Le 
X(1) “19 0 
X(3) Oo 0 


rH OO FF 


O 


o WwW oO x(0) 
goes neem) 
Om Woe) (xe) 
1 o wt | x(3) 
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(ale) 


(2-16) 


(Za) 


(2-18) 


(2-19) 





A signal flew craph of Eq.[ 2-19] is shewn in Figs Vou. 


COMPUTATION ARRAYS 
OA 
Data Array Array 1 Array 2 
Xglk) x4 (k) Xolk) 
x(0) x5 (0) XO) 


<p) 


x9(2) x (L). 





x9 (3) X ( 5 ) 


x9(3) 


miguLeais— “fri signal Flewgraph, N=4. 
Note that tke final results are in bit reversed order with 
respect tc tke disered values er ce Piis-is* silpely che 


scramblir¢e which results frem the FFT algorithno. The 


Ee=Ulct=s are €asily unscrambled by bit reversing, i.ée., 


RACC) X (00) 
a iC) flips or reverses to X (01) (2-20) 
eG 1) X (10) 
meet 1) ery 
Poeeealgcerithy just described is kncwn as the 


Cccley-Tukey Agcrithm or decimation-in-time (DIT) algcritha, 
Since at each stage cf the process the input sequence (i.e., 
time secuence) is divided into smalier sequences for 


processing. 


The kasic building blcck of the algcrithm is the 
"EFutterfly" as shcen in Fig.{ 14]. 
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A X = A+ewk 


8 Y= A-Bwk 
Figure 14 - Decimatior-In-Time Butterfly. 


Each Eutterfly requires one complex multiply (four real 
multipliers) and twe complex adds. The importance cf the 
Evtterfly structure to the FFT flcw graph is that coly ocne 
additicnal memory I1ccaticn is required tc transfcrm an 
N-point sequence stcred in mercry. Thus the intermediate 
results cf tke FFT are stored in the same lccations in which 
the original cata -reing transfcrmed was stored. Since this 
algecrithny stcres toth the input and output sequences in the 


Same locaticn, it is called an in-place algcrithnr. 


As feinted out earlier the results of the DIf algcrithm 
are bit reversed. This means that to get a ncrmal crdered 
O@uteut, tke input data must first be shuffled. Pigs.{ 15a] 
Gtwamh lao.) present the canonic signal flcw gtachs Eor the DIT 
algcritha for the case N=8. Ncte that the data is naturally 


ereemed ir Fig.{j 15a] and is in inverse order in Fig.[ 12bj. 
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Figure 15 - FFI Canonic Ficw Graphs. 





B. CECIFEATICN-~IN-FRECUENCY (DIF) ALGORITHM 


Anctker Forular FFT algcrithm is the Sande-Tukey 
Algcritha or déecimaticn-in-freguency (DIF) algorithn. 


It is 


very sistilar to the DIT algorithm but differs in the 
fcllcwing ways. 


1. 


The input sequence i2i8 is cancer lonee fute 
twe See eee €ach o ength (N/2) samples ir 
tke cClicwing manner. The first sequenc 


ite pe CCNSISts of the first (N/2Z) foints of 
{x(r)}, where as the second sequence es 
cecnsists cf the last (N/2) points cf {x(n)}. 


For a4. given case in the DIT algcrithao the 
input is Fit-reversed while the outfut is ir 
nétural order, whereas the reverse is true for 
the LIF algorithm. In general, however ECtE 
Piceotleand Lie Can gc L£rch nogsal to shuffleé 
cata and vice versa. 


fremecer butterfly -is slightly difterent 
shewr in EietLo de. The difference keing t 
the ccepelex muitipiication takes glace af 
the add-Suktract operation in the 
algcritha. 


as 
hat 
ter 
DIF 


A X= A+B 


k 
Ww 


8 ¥ =(A-8) wk 


Figure 16 - Decimaticn-In-Frequency FPutterfly. 


74 





Taking a lock at the similarities, Ecth algorithms 


require cr tte crder of (N poo Operations tc compute the 
DEI. Ecth algorithms can be done in-place and both need tc 


perfcrm rit reversal at some place during the computation. 


Figs.{15¢cj] and {15d} present the cancnic signal flow 
Craphs fcr the DIF algorithm fer the case N=&. The data is 
nattrally crdered in Fig.{15c] and is in inverse crder in 
mage 15d 5. The two most effective methcds are those 
illustrated in Fiqs.{15b] and {[15c]} since they previde the 
powers cf Wir the correct crder needed fcr computation. 
This elinzinates the need for storage tables. 


The FFI eigorithms examined akove ry nce means exhaust 
aa the rcssibilities. Numerous algorithms have teen 
develcpec, ali with cne basic purpose in mind, and that is 


tc reduce tke ccmputation time. 
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Coeur lLES CSING THE FFT ALGORITHM 


Aas a test cf the FFT algorithm to compute the DFT, CASES 
1,2, and 3, as described in Appendix A, were used as infuts 
for N=3Z. For the CLIT algorithm used, N is reguired to be 
highly ccupcsite and a power cf 2. The value of 32 was used 
so that tke results of this algorithm cculd Ee compared to 
the results cr the CZf and prime transform algorithms. In 
the latter twe algorithms N=31. All figures are drawn with 
ncrralized frequency, but for the sake of discussion, assume 
tbat the input sinuscids have frequency in Hertz anc the 


rectangular cata window "T" is equal te 1 seccnd. 
The funcéemental freguency is 


F = tyT = 1 Hz 
c 


which is tke ninimum recognizable frequency compcnert or 


Eescluticn. the analyzing ctandwidth is 


F = (N/2) (EF ) = 16 Hz 
ma x O 


The output is a one-sided spectra ccnsisting of [ (N/2)+1} 
freguency ccapenents. Figs.{17] and [18] are the results 
fer CASE 1 and those for CASE 2 are shown in Fig.{19] and 
[2C]J. In tke results cf CASE 3, the 22 Hz ccmpfonent aliases 


as a 10 B2 ccafonert. 
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Tit. CHIRP-Z-TRANSFCRM COMPUTER SIMULATION 


A. CZT EROGFAKE DEVELCPEMENT 


A cC@puter Eregram was written to sinulate the C2ZT 
algcritha fcr ccmputing the DFT of an arkitrary wavefcrn. 
This waveicr@ may te real cr cemplex. In carrying cut the 
precessir¢e creraticns, the real and imaginary compcnents are 
represented ry subscrirts "r* and “itt respectively. The 
Simulaticno regram is a result of the follcwing derivation. 


Preceedinac, let 


x-> input sequence 


Keene / 2 
£7 AW te —> ccmplex pre-multiflier (3-1) 
then 
y=xf (3-2) 
yaa eX EK f+ exe. (3=3) 
Eee bp ek a ee 
Sige lifyirg 
faue= = £ - x £ (3=0) 
E oer ed 
vos (ect, te x f) (3-5) 
Zz Eo i 





Now letting 


—-n2/2 
7 7 


* Y 


QO 
tt 
bed 


where "*" denctes convclution 


Seay Avert yYoNgy, + Jy *v + Fy etgV. 
IG JE It 1 x Te 1 Sl 


Ve) 
ul 


Hee) = ey 
Lc i E a a 


Ve) 
i} 


Yoh iv. * jy -*v 
c al Le eae Te 


ioeeemre frcecram this cOnvolution is Situlated 


transversal filter with 2N-1 complex taps v = 
= (N= 


where 


yee eee NS 1) eo. 4 N= 1) 


ana 
W = EXE (-Jj2T/ C) 


Tne third tajcr step in the process, that of 


post-multiplication, is perfcermed Ly letting 


k2/z 
= 4h 


ey 


(3-6) 


(3-7) 


(3-8) 
(3-9) 
(3-10) 


as a 


v 
(N- 1) 


(3-11) 


ccmpiex 


(2=2) 





o. =-'od (3-13) 


ca= Cade + Tq deo+ djqed + j2a a 3-14 
EO 1 a So a ? a i 
Sigsplifyirg 
tea ed. - 9d (3-15) 
i Er la 
mes (ood + gd) (3-16) 
a eae: iS 


Memetne cfecigl case of the DFT, A =1 and & =0 in #9.[1=-23], 
O Cc 


-n 
W=1 and @ =-(1/5) in Eq.{1-24]. Thus A =1, W = exp(-j2 
Cc C 


feeepeeanastc.[s-l], [3-6], and [2-12] can be thought cE as 


ccaplex ezpcnential sequences of linearly increasino 


freguency. They are represented as follows: 


- = ccs (17n2/N) (3-17) 
= +sic (1n2/N) (32-18) 
Ee = cos (tn2/N) (3-19) 
a = sin (mn2/N) | (3-20) 
oe = cos (1k2/M) (3-21) 
oh = -sint (1k2/M) (3-22) 
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Sukstituting [Eq.j3-17] thru [3-22] into EGS. 3-4 ),) (3-5), 
oa ee b=- 10), [2-15], and { 32-16], respectively, describes 
tke imsplemertaticn of the CZT as shown Picterially in 
Pace. ( 23). 


ye 2 pc caims 7) + x a7 8) (3-23) 
Y, = ee 7 8) + BC O = NDE 7.8) (3-24) 
os Ure Coin 737) = Pies one 75) {82729} 
ae = Sane 7 Ny | ee 70) (3-26) 
ae = sce ome) + Drea (3-27) 
i. = Peeks) + 2, Sas aks ee (32-28) 


The sigulaticn pregram is listed in Appendix D ané€ ig a 


morking wedel cf Fig.[ 23]. 
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B. EXAMELES USING THE CZT ALGCRITHM 


momaetect Cf the CZT algorithr to ccempute the DFT, CASES 
1,2, and 3 as described in Appendix A, were used as inputs 
for N=31. Tke imaginary component of the input is set equal 
tc zero since the input consists of only real data. All 
figures ere érawn with normalized frequency, but for the 
sake cf discussicn, assume that the input sinusoids have 
freguency in Hertz and the data windcw "T" is equal tc 1 
seccnd. The data window used for all cases is rectangular. 


Tkhts the furcagental frequency is 


F = i= 1 = 1 = 1 Hz 


= = eee aS oe ee ee 


(ZT) ae (31) (1731) 


which is ir effect the minurum recognizable frequency 
ccaepcnent cr rescliution. fhe analyzing bandwidtk or 


freguency Lrarce is 


E = (NvZ)(F ) = 15.5 #2 
wax O 


However, since the resolution is 1 Hz, the cutput compecnents 
cconsist cf tte D-C ccmponent, the fundamental, and the first 


14 harmcrics fer a frequency range of 15 Hz. 


As in the FFT algorithm, the output of the C2T algcrithnm 
is a cne-sided szpecta consisting cf [(N/2)+1] integer 
frequency ccapenents. Figs.{24] and (25] are the expected 
results fcr CASE 1 and Figs.{26] and [27] are the results 
for CASE z. In the results of CASE 3, it snculd te neted 
ttat the aliasing frequency of 22 Hz fclds back as a J 


ea 





Haypeneccas in the FFT it folded tack as a 10 H2 ccmpcnent. 
This cccurs recause N is odd in the CZT where as in the FFT, 
N is required to be even and highly composite. For the FFT 
arcetithw Ne =22z. 
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IV. PRIME TRANSFORM COMPUTER SIMULATION 


A. ERIME TEANSFORM COMPUTER EFROGRAM DEVELCEEMENT 


So as tc make the mathematics clear, a simple example is 


given tc illustrate the required matrices cf the algcrithn. 


Consider the case N=5 and let W = exp(-j2Z2 7 /N). 
Evaluating the DFT equation directly results ir the 


fcllcwinc N eéquaticnas. 


X(0) = x(0)W° + x(1)w? + x(2)w? + 


X(1) = x(0)we + x(1)wt + x(2)W* + 
a2 = x(0)W? + x(1)Wo + x(2)we + 
mes) = x(0)w° + x(1)w? + x(2)we + 


(14) : 


x(0)w° + x(1)W° + x(2)We + 
Written in matrix form 


X(x) =|F] x(n) 


X(0) wo wl ww 

an fo get yo 
I xc2)} = {wo we we we 
| X(3)| Ww wo wos? 
| x(4)| wO wet we wl 
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x(3)W + x(4)w° 
x(3)W2 + x(4)we 
x(3)we + x(4)w (4-1) 
«(3)W? ~ x(4)wl 
x(3)Wh@4 x(4)wt® 


wo | | x(0) 
wt ly) 
we | | x(2)| 9 (4+2) 
omlex (3) 
wr] | x (4) | 





The matrix {F} is written using the relationships 


nk 


wok . ynk 


mod N 
Wo =] 


For example Ww? = we Since for n=3,k=3 


woe =. W? = exp| (- ae =) (9)| = exp(- j25%) 
(4-3) 
= exp(-J5") = exp [(-j2%)(4)| = wi 
therefore 
X(0) a 1 iL i 1 tne 
X(z)} | 1 0 20owt owe Ow owt | x(1) 
x(2yi={2 we wh wh wl |xc2) (=lp) 
X(3)] a wo who wh ow? | 43) 
x(4)| p20 wh wo wht (4) 
Then deleting X(0) and x(0) q. [bad | is 
X(1) 1 we we wow | x(a) 
X(2) if weoowt owt ow3 Tt | x(a) 7 
= a0.) = 3 i i 2 ( -5) 
X(3) 1! W W W W x(3) 
we ay we wwe wh | x(4) 
x" D F x(n) 
Hor t=5 and r=2, there is a one to one mapping of the 
integers n,k to the integers (n)p,(k)p, i.e., 
eed) 2° = mod 6) 2° = 3 mod 5, 27 = 1 mod 5 
7 
(n)p,(k)p | 2/4)3{1 ae 
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ne 
Matrix {F } is factored into three matrices. 
8 


F = p'cp 


The elements of the {C} matrix are found by using the 
relationship given by Eq. (1-42 |. Thus 


(4-7) 


The elements of the permutation matrix are given by 


Eq.[1-43 |. 


P = © (4-8) 


co 
bt 
© 
O 
oO 


0 QO QO aS 

+ a 0 QO 0 
P = (4-9) 

0 0 aL QO 

0 i 0 0 


A fumetional diagram of the prime transform algorithm 
for computing the Fourier coefficients of real data is as 
shown in Fig.| 30 |. The switch is in the up position for the 
first data sample and is down for the remaining (N-1) sam- 


ples. 


33. 
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B. EXAMELES CSING THE PRIME IRANSFORM ALGCRITHM 


To siwulate the operations of the frime trarsforn 
algcritha, a ccmputer program was written for the odd frime, 


N=31. A listing of the program is given in Appendix E. 


Once again Cases 1,2, and 3 were emplcyed in crder to 
evaluate tte prime transform algorithn. The fundarental 
frequency ard frequency range are the same as that cf the 
G77 algeritbe: for N=31. The results are shcwn in Figs.[{31] 
peu [3c }. AS can pe seen, the results are the same as 


these ccuputed ky the CZT algcriths. 


1C0 





Amen: 


"T eseg ‘WYULTAOSTY WAoJsueTL oWTIg BTA 1OTd epnytusem - 


ZT, OT, a, 9 , 


TE ernst 4 


(pazTTewrzou) aepnytusey 


OT 





"T eseg ‘wWYUITAOSTy waojsuery, eWtid STA 1OTd oseud - ze oeuansty 


Ost- 

nt 06- 
we. 
ae 
9 
Nn 
@ 
I/F if al q ri) = 0 rh wr CH, a 2 fins U a 
ie = I ry cal = ms C3 I = = 4 vs | O O 
0g 
FS 
@ 
U @ 
” 

06 


Ost 


2 





"2 aseg ‘WYYTAOSTY WAOossuerT], sWTIg eTA 4OTd epnytusey - 


S 
3/5 T 2T OT 9 


) 


US 


tS 


fy) 


€€ eansty 


i) 


8° 0 


(peZTTeUr0U) epnyzTusey] 


eS 





EI/ 5 


ci 


HT 


"2 eseg ‘WYUYTIOSTY WIOssUueAL OWTIg BTA 10Td eseud - HE 


[yf 


ZU, On 9 9 t 


SINE Ty 


OSs 


06- 


06 


O8T 


(seerzep) eseudg 


104 





I 


qT 


°€ oseg ‘WUITIOSTy wIOjsueT], swTIg BTA 1OTd opngytusem - GE sansty 


et 0 8, o 4 me 


E L Ze 


9°0 


O°T 


(peZTTeUw10U) apnz use 


105 





5/5 1 


"€ eseg ‘uYyyTAOsTy W1rogsuery, oWTIg BTA OTd eseud - o€ 


ZT oe g _ 9 oN 


INS Ty 


UBis 


06- 


43 0 


06 


OST 


(Seeigep) eseug 


106 





This chéepter addresses the sources cf error ir both 
digital and sampled analog systems with Special attention 


given tc the sampled analcg implementaticn of tte C2T 
algceritha. 


A. CIGIIAL FEIT SYSTEM 
In a digital FFT system, the sources cf errcr tay be 


Suagparizedc as fcllews: 


igeeGUiantization Of the data at the input A/E 
ccrverter. 


Ze Inaccurate transformation due te errors an 
firite representaticn cf the coefficients os 
inetoue butterflies. 


BeaeeSeLodoti nolse in truncating the result of 2a 
gultiplication and scaling the data tc 
prevent cverflow. 


An «ms guantization error in the AYE ccnverter can be 


defined as 
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where £1 is the nugber of bits used to represent the Poipu 
cata. This ncise dces not scale with signal size and 
dcginates cnly at low signal levlels. 


The errer due to the inexact representaticn cf the 


yen: ee k 
multiplier ccefficients ae Can be expressed as the ratio of 


mean square cutput error to gean square cutput sigral to 


give 


~ (2) (3 
B. = [reg,s\2 sda (5-2) 


where bi ais the number of bits used to represent the FFT 
ccefficients. This eguation predicts that the variance of 
tke errcr increases at a rather slow rate with N. 


A measure of the rms noise cutput to the rms signal 
output due te errcrs produced by rcundcff and signal 
trtncaticn te prevent overflcws in the butterflies is given 
by 


73 
Tus(@rrcr) \/n2 (0.3) 8 (6-3) 
rms (signal) rms (input) 


where £2 is the nugker of bits used to represent the result 
Sumetavce cperation. This upper bound of the errer~to=signal 


is seen tc ircreasée asvV/N, and is the mcst important scurce 


GmeerEcr if @ digital FFT. 
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EB. SAMPIEL ANALOG CZT SYSTEM 


Instead cf errers resulting from registers containing a 
finite werd length, as in the digital FFT system, the cutput 
accuracy cf the sampled analog CZT system is limited by the 
cCl's and SA&k's Sensitivity to the enviroment (temperature, 
hubidity, etc.,) and fabricaticn imperfecticns. The scurces 
cr errer cf a CCD CZT system using pfre~ and pfest-multirlying 
ccefficients stored in either a ROM or PROM are as fcllows 


T Thergeal noise-> this, source cf error is 
apelegous to Quantization in a digital FFI 
recatuse it ‘generates an errer “which is 
independent of signal level. 


2. WTIke errers introduced by inaccuracies in the 
pre-and fost-multipliers. 


3. Ike errors introduced by the inaccuracies in 
tte tap weights of the transversal filter. 


4, Ike error due to mismatch and nonlinearities 
in tte differential current amplifiers. 


5. Tke error due to mismatches of the 
pcst-multiplier outputs. 


The exrzrers addressed in this study where those due to 
imaccuracies in the pre- and post-multipliers and the 
inaccuracies in the tap weights cf the transversal filter. 
The errors ir the pre- and pcst-multipliers were modeled as 
being a percentage "p* of the theoretical ccefficient value 


m(D), 1-€., 


We EL m(n) J | (5-4) 
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and to te uriformly distrituted on (-%/2, %/2). The 
transversal filter was modeled as a tapped analog delay line 
with the tap weights determined by the value of the external 
Mesistors. Ihe eé¢rror in the value of the tap weight was 
mcceled as keing a percentage "p" of the theoretical tarp 
eriome value E(n), i.e., 


B = p{h(x),J (5-5) 


ard tc be urifcrmly distributed cn (- 872, B/z). These error 
mcdels were included in the simulation progran [Appendix Dj}. 


See oeowelLIPWITY TESTS 


Three cifferert input wavefcrms sere used to evaluate 
the perfcrmace of the CZT algorithm for the fcllowing 


percentaces 
p= C#, 4%, 10%, 20%, 30%, and 40% 


In all three tests N=512 and a rectangular data windcw equal 
tc 4 seccnds was used. The value of 4 seccnds was used for 
ease cf flctting the output spectrum. Thus the freguency 


rescluticn is fff = 0.25 and the frequency range is fyf = 
iS iS 


€4.C. 


The first waveform consisted cf the impulse furction 
shewn in Fig.{37]. The time impulse responses of the ccsine 
anéd sine filters of the convolver for _ the different 
percentages are shown in Figs.) 35 jerome oon. FOr 
p=(%,Fig.{£11 shows the expected results, i.e., egual 
macnitude frequency components covering the entire frequency 


rance. As tte percentage of error increases, examination of 
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the figtres reveals that the average vale cf the magritude 
of the frequency ccmgonents is greater than 1.0. 
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Fig.{57] was tke second waveform used fcr computin¢ the 
DEI. TkEe magnitude of the one-sided Spectra 1s shcwn in 
Figs.[58] thru [€3]. Degradaticn of the Outfut spectrum was 
hever sc severe that the (sin x)yx function was not 
reccgnizatle, tut as pointed out above, the average value of 
the magnitude shows a moderate increase as the value cf "pt 
increases. 


The third infgut waveform, shown in Fig.{64], ccnsisted 
of the fcllicwirg ccmgenents. 


Amf. ae Phase 

Saa ft 

sine tare 14.125 Qo 
SGS1£e a ae 15.0 Qo 
Sine ieee 16.0 Qo 
cosine Uae 16.25 Qo 
Sine -0.4 1735 60° 


The ncrmalized one-sided spectra and fhase plete are 
eucmio tr Figq.{65} thru {76}. In Fig.[65], the spreading of 
the 14.125 ccmponent is due to leakage. the 16.0 and 16.25 
ccuponents e@pprear to be indistinguishable, hcwever, this is 
a iss-leading result caused Ly the continucus curve flot. 
If a point ocrafh had been used, the two ccmfonents would be 
readily identified, an expected result since the freguency 


rescluticr is fff = 0.25. 
iS 


Exaginaticn of the phase flcts' reveals that the 
degradaticn is mcre severe than that of the fagnitude flcts, 
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even for sgall values of "p". As in the previous exaugles, 

the magnituce flcts also exhitit a general increase in the 

cCmpcnent magnitude except for the f/f = 16.0 compcnent. 
Sc 


= 


{his ccufgenent <shows a very distinguishing decrease as "tc" 


increases. Ey re-examining the magnitude flots of the 
impulse waveform, it is obvious that the f/f = 16.0 
s 


ccmponent will show a decrease as "pf" increases. These 


Flcts alsc show that a more severe degradaticn cf the 
cne-sided spectra would result for a similar set Gt 


ccopcenents near f/f = 48.0. 
Ss 
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VI. CONCLUSIONS AND RECCNMENCATIONS 


The FFT ¢lgcrithm is a well established method for 
ccuputing the DF1. The chirp-z-transform algorithm can also 
ccmpute the DFT and since the bulk cf the recuired 
cc@jrutaticns can be performed by a transversal filter, it 
lends itself naturally to implementation with sampled analog 
devices. Ike prigwe transform algoritimhas demonstrated its 
akility tc ccmpute the DFT and like the CZT algorithm, can 
tke iaplemerted with sampled analog devices. The fact that 
the multipliers required by the CZT algorithr are refglaced 
by permutinc memcries in the prime transfcrm algorithn may 
cr may nct re an advantage. It will deperd upeon the speed 
and dynamic range required for a particular applicatiocn. 


This paper evaluated the errors in the weights cf a 
transversal filter for weighting provided externaliy tc the 
CCL ty resistors. Preliminary results indicate that the tap 
weights may vary as much as 20% before the frequency 
ccmpcnents are greatly distcrted. HOWeVEr, additicnal 
aralysis neeés to te conducted before conclusive results are 
presented. In additicn to resistor weighting, a Split-sate 
weight technique has teen developed which weights the taps 
of the transversal filter during fabrication. For this type 
cf weighting, quantization errors are introduced. This 
suggests that a different errer model will re required to 


cceréuct a sersitivity analysis 
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APPENDIX A 


ALGORITHM TEST DATA 


Several test waveforms were computer generated to 
evaluate the ferformance of the varicus alacrithms used to 
ccupute tke [FT. 


Case 1-> This waveform ccnsisted cf the fcllowing 
ccrfpcnents: 


AMP. f/t | Phase 

Shift 
Ea 0.8 0.0 Qo 
cosire feo 2.9 go 
Sire 2.0 6.0 go 
ccsire 1.0 14.0 450 


The intert cf this case was to see how well the algorithms 
identifieé tke sagnitude and phase infcrmation cf each 


frequency ccayrcnent. 





Casé-> This waveform ccnsisted OL the fellowing 
ccipenents 


Amp. a Phase 

Shift 
rao Gee 0.0 Qo 
cosine 1.¢ AL go 
Sine 2.¢ 6.0 Qo 
sire el etc UO 5 Qo 
ccsire aren 14.0 450 


This case illustrates what happens when leakage is present. 
The data nindcw truncates the fyf = 10.5 component at a 
s 


non-integer s#ultiple of the wavelength, which is the cause 
cf leakace. 


Case 3-> This waveform ccnsisted of the fclicwing 


Gceijecnents 


Amp. E/E Phase 
s 

Shift 
c-C 0.8 0.0 Qo 
ccsine 1.€ 2.0 Qo 
sire 2.6 6.0 Qo 
ccsire 1.0 14.0 459 
ccsire -2.C 2250 Qo 





This wavefcrm has a component greater than (f )/2, thus 
s 
this case was used te illustrate the aliasing phenorencn and 


the effect N (being odd or even) had on the folding cf this 
ccmpenent inte the analyzing Ltandwidth. 
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APEFENDIX B 


ERIEF REVIEW OF INTEGERS AND INCICES 


The fclicwing definition and theorems are provided fcr a 
better understanding of the derivaticn of the prime 
transfers alccriths. 


B VDEFINITICNS 


Unit-> a unit is an integer that divides every integer. 
Since +1 and -1 divide every integer, they are both units. 


Null elewent-> it is an integer that divides cnly 
itself. Nc integer different frem zerc is a null element. 


Zero is the nuli element of the rational integers. 


Associates of an integer-> they are the results of 
multiplyirg the integer by the units. 


Prime-> it is an integer, not a unit, that is divisible 
by cnly ite asscciates and the units. This definition 
implies that the greatest common diviscr of a prime "p" and 
an integer "a" is 1, or the pesitive associate of "p". 


Compcsite-> it is an integer that is not the aull 
element, a trit, cr a prime, e.g. the integer 256 is said to 
be highly ccapesite while the integer 563 is moderately 


ccmpcsite. 
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Relatively prime-> two or more integers are frire to 
each ctker cr relatively prise, if their greatest ccumcn 


divisor is +1. The integers 6,-9,14 are relatively prime. 


Ictient cf m-> denoted g(m) is the number of positive 
integers less than "m" and relatively prime to "nm". 


EB. THEOKEES 


ie Iwo integers "a" and "b" are said to be ccncruent 
modulo ard integer "mn" if "m" is an integer divisor cf a-b. 
This relaticr btetween "a" and "b" is mest often denoted Ly 


azb mod (m) | (E-1) 


ané is reaé "a is congruent to b mod (m)" Eut for the ease 


cf ncetaticn, it is denoted as 
((a4)) = b mod (m) (E>2) 
in the derivation cf the prime transform algcritha. 


If "at is congruent to "b" modulo (m) and OSbsa,"b" is 
called the residue of a modulo (m}. TE “m" ais net an 
integer diviscr of a-b, "a" and “b" are said to be distinct 


mpcaulo om. 


Ze If "at and "m"™ (<0) are relatively prime integers, 


then 


Mae = 1 mod (am) (E-3) 


158 





Corollary of 2. If"p" is a prime which is not a 
divisor cf ar integer "a", then 


ny 


= 
((@)) = 1 mod (p) (E~4) 


If "a" and "a"(>0) are relatively prime integers, and "u" is 
the least fositive integer for which 


u 
((@ )) = 1 mod (u) (E-5) 


“u" is called the expenent to which "a" belcngs modulo (m). 
If the expenent to which "a" belcngs mcdulc (m) is g(m), “a" 
is definec as a primitive roct modulo (nm) (Or "a ~ pearl tive 
eect CL eg). 
3. The cnly integers which fossess primitve roots are 
D E 
Pewee gee , “Shere "pp" is an cdd prime. if “m" does have a 
Erimitive rect, it has g(g(m)) distinct rocts module (nm). 
Let "p" ke any prime, and let "r" be any primitive root 


cf£ "p", To each integer "a" relatively prime to "rf" there 


ccrrespends a unigue integer "i" such that 


i . 
((a)) = xc mod (p) (0<i<p-1) (E-6) 


Mbe integer “i is called the index to the base “rN cf a 


fedulo (f). This is written as 


i = ind a (E~-7) 
rE 
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Prcperties of Indices 


{a) ind 1 = 0 


: (E-8) 
() ad tent = rt Dia die {B-9) 
(c) ind (ab) = ind a + ind b mod (p-1) (B-16) 
c rc c 
= ° Q « 
(dq) ind a@ = n ind a mod {p-1) {B-11) 
r c 


The power residue of an integer “b" to the base "xr" is 
Simply the integer "s" such that 


b 
({(s)) = © mcd (p) (0<s<p) (B-12) 
Since tke power residue of inda is ccngruent tc a 
c 


mcdulc (p), the concept of indicies and fewer residues is 


analcgous te that cf logarithms and antilogrithms. 
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APPENDIX C 
OF THE FFT ALGORITHM 


PROGRAM LISTING 
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